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■ i—{ , 1 INTRODUCTION 



ABSTRACT 

We investigate gravitational lensing by a realistic model of disk galaxies. Most of the 
mass is contained in a large spherical isothermal dark matter halo, but the potential 
is modified significantly in the core by a gravitationally-dominant exponential disk. 
The method used is adapted from a very general multi-lens ray-tracing technique 
developed by Moller. We investigate the effects of the disk-to-halo mass ratio, the 
disk scale length, the disk inclination to the line of sight and the lens redshift on 
two strong-lensing cross sections: the cross section for multiple imaging, and the cross 
section for large magnifications, in excess of a factor of 10. We find that the multiple- 
imaging cross section can be enhanced significantly by an almost edge-on Milky Way 
disk as compared with a singular isothermal sphere (SIS) in individual cases; however, 
when averaged over all disk inclinations, the cross section is only increased by about 
fifty per cent. These results are consistent with other recent work. The presence of a 
disk, however, increases the inclination-averaged high-magnification cross section by 
an order of magnitude as compared with a SIS. This result has important implications 
for magnification bias in future lens surveys, particularly those in the submillimetre 
waveband, where dust extinction in the lensing galaxy has no effect on the brightness 
of the images. 

Key words: methods: numerical - galaxies: fundamental parameters - galaxies: 
spiral galaxies - cosmology: theory - gravitational lensing 



X 



There is ever increasing interest in gravitational lensing by 
individual galaxies. Deep surveys in the radio waveband, 
for example the CLASS survey (Myers et al. 1995), have 
identified many galaxy lensed systems. In many cases, mod- 
els of the lensing mass distribution as spherical or ellipti- 
cal isothermal halos have been reasonably successful in re- 
producing the observed image positions and magnification 
ratios (for example, Rhee 1991). An increasing number of 
observed lens systems have, however, image configurations 
that are not consistent with this simple picture (for exam- 
ple Jaunsen & Hjorth 1997). These have recently lead to 
the consideration of lens models that include a disk compo- 
nent within a dark halo (Bartelmann & Loeb 1998; Keeton 
& Kochanek 1998; Koopmans, de Bruin & Jackson 1998; 
Mailer, Flores & Primack 1997; Wang & Turner 1997). 
Since a disk with mass comparable to that of the halo would 
be gravitationally unstable, the mass contained in the disk 
is not expected to exceed more than about 10 per cent of the 
halo mass. Even so, within a few kpc of the core of such a 
galaxy the projected surface mass density on the lens plane 
is dominated by the disk component and this enhancement 
should lead to an increase in the strong-lensing cross sec- 



tion. The effect is expected to be particularly significant for 
a nearly edge-on disk. 

Wang & Turner (1997) calculated the cross sections for mul- 
tiple imaging for lenses consisting of a singular isothermal 
sphere (SIS) and a thin uniform disk, and found a large in- 
crease at large disk inclinations; however, the model of a 
uniform galaxy disk is not very realistic. Mailer et al. (1997) 
used an exponential disk mass distribution to fit the two- 
image lens system B 1600+434 (Jackson et al. 1995) with 
reasonable success. Keeton & Kochanek (1998) presented 
a model that combined three different elliptical mass com- 
ponents to produce an analytical approximation to an ex- 
ponential disk potential with a flat rotation curve at large 
distances. This model involves many parameters and the 
potential differs significantly from that of an exponential 
profile near the core of the galaxy. While we were preparing 
this paper, Bartelmann & Loeb (1998) discussed the statis- 
tical significance of lensing by spiral disks using a technique 
similar to ours to derive strong lensing optical depths, in- 
cluding the effect of dust extinction along the lines of sight 
to different images. We consider these effects in detail in a 
subsequent paper (Blain, Moller & Mailer; in preparation). 

This work was motivated by interest in the effects of 
disk lensing on the population of lenses in the millime- 
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Figure 1. Rotation curves for an exponential disk in an SIS halo 
with ov = 150kms — 1 . The circular velocity v c = V2a v . 



tre/submillimetre waveband, where, in some circumstances, 
magnification biases are expected to be very large and up 
to 5 to 10 per cent of sources detected in a survey could be 
lensed (Blain 1996, 1997, 1998). An enhanced strong-lensing 
cross section due to disks could increase this effect even fur- 
ther (Blain et al. in preparation), especially because dust 
extinction is unimportant in the submillimetre waveband. 

We adapted an existing general ray-tracing method, 
which was first developed to study weak lensing (Moller 
1997), to solve the problem. The method and lens model 
are outlined in Section 2. We then investigate the effects of 
the surface density and scale length of the disk on the cross 
section for the formation of multiple images, which is es- 
sential for determining whether a source in a survey will 
be classified as a lens in Section 3, and the cross section 
for large magnifications, which determines the importance 
of magnification bias (Borgeest, Linde & Refsdal 1991) in 
Section 4. In Section 5 we compare our results with those 
of previous authors and briefly discuss their implications 
for lensing statistics and magnification biases. We assume 
that the density parameter Oq = 1 and Hubble's constant 
-ffo = GOkms" 1 Mpc -1 . Angular diameter distances are cal- 
culated assuming a smooth distribution of matter. 



2 METHOD 

Even for a uniform disk, the analytical derivation of the form 
of caustics and critical lines is very involved. In order to 
study lensing by exponential disks, we adapted an existing 
general multi-lens-plane ray-tracing routine (Moller 1997), 
in which a regular grid of Ri x Ri images is mapped onto 
the source plane by propagating light rays through a set of 
N lens planes. In the present computations, N — 1. The 
simulated source plane is finite, and so care was taken that 
no rays are deflected beyond its boundaries. In principle, 
extended sources can be treated, but all the results are de- 
rived here for point sources; the brightest regions observed 
in far-infrared-luminous galaxies at low redshifts are typi- 
cally less than a few hundred parsecs across (Solomon et al. 
1997), corresponding to an angular size of order 10 -2 arcsec 
at redshift z — 2, and so this approximation is unlikely to 



affect our applications significantly. To generate maps of the 
total magnification on the source plane, we find all the im- 
ages of a regular grid of Rs x Rs point sources using the 
'grid search method' described by Schneider, Ehlers & Falco 
(1992). Once all the images of each source are located, we 
find its total magnification by summing all the individual 
magnifications. 

The accuracy of this procedure depends on the resolu- 
tion of the grid on the source and image planes. If the image 
resolution is too coarse, then the grid search method may 
miss images near caustics, and so the total magnification is 
underestimated. This problem can be partially overcome by 
making Ri 3> -Rs • There is a trade-off between the accuracy 
of magnifications near caustics, the coarseness of the result- 
ing map and the computing time. The cross sections for dif- 
ferent magnifications are determined by adding together the 
areas of all pixels in the source plane that contain sources 
with that magnification. Hence, a small value of Rs leads to 
an overestimate of the high-magnification cross section as 
the average magnification over the whole pixel area is not 
sampled adequately. To verify the accuracy of our proce- 
dure and to find optimum values for Ri and Rs, we tested 
its predictions against analytic results derived for a SIS lens. 
For resolution ratios Ri/Rs > 2 the results agreed with the 
analytical predictions to within about 5 per cent. In all com- 
putations, Ri = 1000 and R s = 400. 

The deflection angle a is calculated by adding together 
the deflection due to the two different lens components, a 
SIS halo and a disk, 



a = otsis + ctdisk- 



The deflection due to the halo, 



4GM 
asis = p 2 r, 



(1) 



(2) 



where Mo is the mass contained within radius Ro and r is 
the impact parameter. The projection of an inclined thin 
disk onto the lens plane produces an elliptical surface mass 
density profile. We evaluate c*disk numerically using a useful 
elliptical surface potential formalism developed by Schramm 
(1990). 

We consider a single model lens at redshift z\. The 
lens consists of a SIS dark matter halo with velocity dis- 
persion <7 V and a disk with a surface mass profile E(r) = 
Eo exp (— r/r s ), in which r s is the radial scale length. For 
clarity, no bulge component is included, but our model can 
incorporate additional mass components without difficulty. 
The effect of including a bulge is discussed briefly in Sec- 
tion 4. In Fig. [j] four realistic fiat rotation curves predicted 
by the mass models used in this paper are shown. Values 
of cr v = 150 km s" 1 , E = 1O 9 M kpc" 2 and r s = 3.2 kpc 
produce a reasonable fit to the rotation curve of the Milky 
Way (Binney & Tremaine 1987). 



3 MULTIPLE-IMAGING CROSS SECTION 

For a SIS lens, the cross section for multiple imaging, 

2 



sis 



S^D L s 



(3) 
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Figure 2. Multiple-imaging cross sections for a lens consisting of a SIS halo, with a velocity dispersion ov = 150kms — , and an 
exponential disk with a central surface density So and radial scale length r B . In (a) So is fixed at 10 9 Mq kpc -2 : the upper, middle 
and lower curves correspond to r B = 8.0, 5.0 and 3.2 kpc respectively. In (b) r B is fixed at 5.0 kpc; the upper, middle and lower 
curves correspond to So = 2.6 X 1O 9 M0 kpc -2 , 1.0 X lO 9 M0kpc -2 and 0.4 X 1O 9 M0 kpc -2 respectively. In all cases z\ = 0.2, 
z s = 1.5. 



(Schneider et al. 1992); Dls is the lens-source angular diam- 
eter distance. If the lens redshift z\ = 0.2, the source redshift 
z s = 1.5 and cr v = 150 km s -1 then a^ s — 15.1 kpc 2 . Includ- 
ing a disk component increases <r m if the inclination angle of 
the disk 9 exceeds a critical angle 9 C (equation 10 in Mailer 
et al. 1997), 9 = and 90deg corresponding to a face-on 
and edge-on disk respectively. Note that we add the disk to 
a halo with a fixed velocity dispersion and consequently, our 
pure SIS model does not take into account the disk mass. 
We use this approach because we aim to demonstrate how a 
disk can increase the lensing cross section without requiring 
a more massive halo. For a Milky Way type galaxy, the disk- 
to-halo mass ratio within a ten kiloparsec radius is about 1:2, 
and including the disk mass into a pure SIS model would in- 
crease the velocity dispersion of the halo above the observed 
value. 

The variation of <r m as a function of inclination is shown 
in Fig. |2| for various disk parameters. In all cases, <r m in- 
creases monotonically with increasing inclination, and the 
effect is considerably more significant for more massive disks. 
Although the parameters are chosen so that the upper curves 
in Figs.^(a) & (b) are derived for disks of the same to- 
tal mass, the predicted cross sections differ significantly at 
small inclinations: the disk with the larger surface mass den- 
sity has the more significant effect. These results agree with 
the expected dependence of 9 C on the surface mass density 
reported by Mailer et al. (1997). The effect on the cross sec- 
tion ratio of changing z\ is negligible, as compared with the 
dependence on the surface mass density and scale length. 

The multiple-imaging cross section is increased sig- 
nificantly at large inclinations, but if r B = 3.2 kpc and 
Eo = lO 9 M0kpc -2 , the inclination-averaged multiple- 
lensing cross section, 



tr/2 



i9am(9)d9, 



(4) 



only exceeds that of the SIS halo by about 50 per cent. 
Hence, although high-inclination disks can increase the 
multiple-lensing cross section significantly as compared with 



a SIS, the inclination-averaged effect of disks is small: the ef- 
fect of a disk can be very important when modelling the im- 
age configurations in individual lens systems, but the pres- 
ence of a population of disks should not have a significant 
effect on the number of lenses that will be identified in a 
large sample of candidate sources. We also investigated the 
relative number of images. For a Milky Way type galaxy we 
find that the inclination averaged cross section for the for- 
mation of double and triple images is approximately equal, 
with about 40-50 per cent of the total multiple image cross 
section being provided by each. Only about 10 per cent of 
multiply imaged sources will be split into four images. The 
cross section for a source being lensed into five images is 
zero within the numerical error. Note that these results do 
not take into account magnification bias. 



4 HIGH-MAGNIFICATION CROSS SECTION 

For a SIS lens, the cross section for magnifications fi > A, 



_SIS/ 



>(A-1) 



for A > 2; 
_2 , forl<A<2, 



(5) 



as can easily be obtained from the magnification curve 
(Schneider et al. 1992). For large magnifications, however, 
we predict that <r M is much larger for a SIS plus disk lens 
as compared with a SIS alone. Magnification maps derived 
at inclinations of 60, 70 and 85 deg for a lens similar to the 
Milky Way are presented in Fig.^, and clearly show the dis- 
tinctive shape of the caustic, with an elliptical component 
due to the halo and a diamond shaped component, the "as- 
troid", due to the disk. For large inclinations, the size of the 
astroid increases significantly, leading to the increase of the 
multiple-imaging cross section cr m discussed above. 

The cross section ratio a^/a^ for magnifications 
greater than A for a galaxy lens similar to the Milky Way are 
presented in Fig. 4. The cross sections are derived by adding 
up the area of all pixels on the source plane with magnifi- 
cations greater than A. The calculated ratios become less 
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Figure 3. Magnification maps for a lens centred at the origin that consists of an SIS halo with <r v = 150 km s -1 and a disk with 
r s = 3.2 kpc with its major axis lying along the x-axis; z\ = 0.2 and z s = 1.5. The grey-scale represents the total magnification of 
a source located at x,y in the source plane. In (a), (c) and (e) So = 10 9 Mq kpc - 2 . In (b), (d) and (f) So = 3 X 10 9 Mq kpc~ 2 . 
6 = 60deg in (a) and (b), 9 = 70dcg in (c) and (d) and 9 = 85deg in (e) and (f). The outlined box in the bottom-right corner of 
each panel encloses a constant area and is shown to indicate the different scales. 
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Figure 4. High-magnification cross sections for a galaxy with a 
SIS with <Tv = 150kms — 1 and a Milky Way disk. To indicate the 
uncertainties involved the curves are not smoothed. 

smooth at large magnifications because the number of high- 
magnification pixels on the source plane is relatively small; 
only about 20 pixels in the source plane have magnifications 
larger than A = 80. The magnitude of the fluctuations in 
the ratios can be used as an estimate of the uncertainty in 
the results of these computations. A large increase in the 
cross section ratio is predicted for inclinations that exceed 
60 deg, at magnifications greater than a certain threshold. 
This threshold magnification is a strong function of inclina- 
tion, and is reduced as the inclination increases. The cross 
section ratios shown in Figs. 4 & 5 also have a characteristic 
peak beyond this threshold magnification. Hence, although 
the largest cross section ratios are predicted for almost edge- 
on disks at magnifications greater than 50, at magnifications 
of about 10 a disk at 9 = 70 deg is expected to produce 
a larger cross section than a disk at 90 deg. Hence, nearly 
edge-on disks are not necessarily more effective in produc- 
ing large magnifications as compared with disks at smaller 
inclinations. This is because points interior to the astroid 
caustic lie closer to caustic lines and experience larger av- 
erage magnifications when the caustic is smaller, as it is at 
moderate inclinations. This effect is clearly demonstrated 
in Fig. 3; although the area of the astroid caustic is larger 
in the lower panels, the mean magnification in the enclosed 
region is smaller. 

The inclination-averaged high-magnification cross sec- 
tion <t m shown in Fig.^ was derived by evaluating <r M at 
five-degree intervals from to 90 degrees and approximating 
the integral (equation 4) by a sum. The result is very signif- 
icant: <J M is greatly increased by a randomly-oriented galaxy 
similar to the Milky Way as compared with a SIS lens, and 
o'/csis ~ 10 at magnifications greater than about 50. This 
increase could have a very significant effect on magnification 
bias in galaxy surveys. 

The dependence of <r M on the scale length and surface 
density of the disk is shown in Fig. 5. As in the case of 
the multiple-imaging cross section, the more massive disks 
produce a more significant effect. Again, if the central sur- 
face mass density is increased, then the relative significance 
of disks with smaller inclinations increases. Adding a cen- 
tral bulge does not have a significant effect on the high- 



magnification cross section; for a spherical bulge with a r 1//4 - 
law density profile, extending out to one kiloparsec from 
the centre and containing a total of 10 per cent of the disk 
mass, we find that the high-magnification cross section is in- 
creased by no more than 50 per cent. The projected surface 
mass density profile for an elliptical halo is similar to that 
for an inclined disk. Thus, we expect that a non-spherical 
halo which is flattened along the disk is likely to enhance 
the lensing effect of an inclined disk. However, this effect is 
not expected to be significant. 

The high-magnification cross section does not depend 
strongly on the redshift of the lens. This suggests that the 
high-magnification cross section in both, the disk plus SIS 
halo model and the pure SIS halo model scales with the 
angular diamater distance in a similar way. Thus, we do not 
expect o",j/<t^ is 1° be strongly dependent on the value of the 
density parameter Qq or the comological constant A. 



5 DISCUSSION 

We have studied the effects of lensing by an exponential 
disk inside a SIS on two important parameters; the multiple- 
imaging cross section cr m , and the high- magnification cross- 
section a fl (A). The lens model we use has a density profile 
similar to that of the Milky Way. Both <r m and ff^(A) axe 
sensitive to the lens parameters, in particular to the inclina- 
tion angle and the disk-to-halo mass ratio. 

The results indicate that disks, and especially almost 
edge-on disks, have a significant effect on the image configu- 
rations produced by individual lenses, which agrees with the 
results of previous work by Mailer et al. (1997) and Wang & 
Turner (1997), but that the statistical effect on <r m averaged 
over all disk inclinations is unlikely to be very important. 

The more significant result is that including an expo- 
nential disk at an inclination greater than about 60 deg in- 
creases cr M at A > 10 by an order of magnitude as compared 
with an SIS lens, and that the inclination-averaged cross sec- 
tion <r M is increased by a similar factor. Hence, disks could 
have a large effect on the predicted abundances of strongly 
magnified galaxy lenses. This implies that in surveys for 
which dust extinction is unimportant, such as those in the 
radio and submillimetre wavebands, the number of detected, 
strongly lensed sources will increase (Blain 1996, 1997), and 
that they will tend to be preferentially associated with edge- 
on spiral galaxies (Blain et al., in preparation). The prelimi- 
nary results of the CLASS survey suggest that this is indeed 
the case (Browne et al. 1997). A recent preprint by Bartel- 
mann & Loeb (1998) discussed the statistical significance 
of lensing by spiral galaxies using a technique similar to 
ours and included the effect of dust extinction in the lens- 
ing galaxy. They did not discuss the effects of different scale 
lengths or disk masses on the high-magnification cross sec- 
tions and so it is difficult to compare our results in detail, 
but they appear to be in broad agreement, bearing in mind 
that they assume Q.q = 0.3 and Ho — 50kms _1 Mpc^ 1 and 
that their disk model deviates from our exponential profile. 



6 CONCLUSIONS 
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Figure 5. High-magnification cross sections for galaxy models with different disk parameters from those assumed in Fig. 4. In (a) 
and (b) the mass of the disk is 2.5 and 3 times larger as compared with the model in Fig. 4. 



(i) The configuration of images in an individual lens can 
be modified significantly by the presence of a disk compo- 
nent at a large inclination, but when averaged over all in- 
clinations the multiple-imaging cross section of the lens is 
increased by only about 50 per cent. 

(ii) The large-magnification cross section is increased con- 
siderably by including a disk component in the lens model. 
For a Milky Way disk the effect is very significant at incli- 
nations in excess of a moderate value of about 60deg, and 
so the inclination-averaged cross section is also expected to 
increase significantly, by an order of magnitude at magnifi- 
cations of order 10 as compared with a singular isothermal 
sphere lens. 
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